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We investigate an effect that an ideal Lorenz number Li of correlated metal shows peculiar Fermi 
surface dependence, which is caused by the onset of a particular channel of Umklapp scattering. We 
evaluate Li for some simple models and transition metals, and note that Li for Na^CoC^ decreases 
sensitively as x approaches an Umklapp threshold around x c ~ 0.6. 
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I. INTRODUCTION 



The discovery of large thermopower in Na^OoC^ [H has 
prompted quest for related good thermoelectric (TE) ma- 
terials and encouraged experimental and theoretical re- 
search on TE properties of strongly correlated electron 
systems. In practical TE application, the important 
material parameter is the dimensionless figure of merit 
ZT = S 2 /(k/(jT), where S is the Seebeck coefficient, n 
the thermal conductivity, a the electrical conductivity, 
and T the absolute temperature. Hence it is equally as 
important to enhance S as it is to decrease the Lorenz 
number L defined by L = k/ctT^. The latter however 
is usually hampered by the Wiedemann-Franz (WF) law, 
according to which L should be a universal constant, e.g., 
L = ^(k/e) 2 for impure metals. In fact, normal metals, 
which generally have low S, will not be good TE devices 
with ZT > 1, unless the WF law is overcome. It would 
thus be interesting to investigate a possible material de- 
pendence of L to see if it is a controllable variable in 
principle. 

We discuss the ideal value Li (which is simply de- 
noted as L below) of a correlated electron system on 
a rigid lattice without impurity in the low temperature 
limit for the purpose of elucidating its material spe- 
cific dependence. Indeed, for pure transition metals, the 
Lorenz numbers have been observed to vary from metal 
to metal 0,13, IE Theoretically, Herring had derived 

a universal constant Q, while the problem had also been 
addressed specifically by treating s-d hybridization as 
perturbation, namely, by a two band model of conductive 
s electrons scattering off localized d states [E flOj - Her- 
ring's argument based on the Fermi liquid theory should 
be the proper approach to the problem at low temper- 
atures. We derive a dimensionless factor (Eq. (fT0|) ) so 
as to modify his constant result, with which we investi- 
gate the Fermi surface (FS) dependence of L for some 
typical cases concretely. After showing nontrivial results 
obtained for simple models, we give results of numerical 
evaluation of the Lorenz numbers of some transition met- 
als. Lastly, we find it interesting to investigate a model of 
the sodium cobalt oxide Na^CoC^, because a nontrivial 
x dependence of L is expected theoretically owing to its 
simple cylindrical FSjllj. 



II. FERMI LIQUID THEORY 

We begin with the linearized transport equation of a 
Fermi liquid under the temperature gradient VT and the 
electric field E, 



T 



VT • v„ + eE • v, 



m, (i) 



TI[<p] = 



p f ,k 



'Wi 



VP 



i p n p ,(l - n p _ fc )(l - n p , +k )(ip p + ip p > 



~<Pp-k 
where n p = n{e p ) 



f P '+k)S(e p + 

■ l/( e (e P -At)/T. 



1) is the Fermi distri- 



bution function, and Wj? p represents the transition prob- 
ability of quasiparticle scattering; p,p' — > p — k,p' + k. In 
terms of the solution (p p of Eq. (JTJ), the electric and heat 
currents carried by quasiparticles are respectively given 
by 



J = 2eY^v p - — ip p , Q = 2 



P 



Vp(£ p - fi) — ip p , 



where e p and v p are energy and velocity of quasiparticles, 
respectively. In comparison with the phenomenological 
formulae J = ctE -ctSVT and Q = TctSE-koVT, we ob- 
tain the transport coefficients a, S, and kq, the thermal 
conductivity at zero electric field. The thermal conduc- 
tivity at zero current, k, is given by k = kq — TaS 2 ^. 
Formally, the above results expressed in terms of the 
rcnormalizcd quantities e p and v p bear resemblance to 
those of a weakly interacting Fermi gas. Nevertheless, 
many body effects are included not only in the renor- 
malization for the individual quasiparticle, but also in 
the field induced shift of the quasiparticle distribution 
tp p , which implicitly includes a collective Fermi liquid ef- 
fect depending on Landau parameters |l2l. If3| . We do not 
write down the explicit expressions for them as they are 
irrelevant for our purposes in what follows. 

For definiteness, let us assume that the currents J and 
Q flow in the x-direction. Then, for the two functions h 
and ?2 defined by ip p = —eE x l\ + d x Tl-z, we obtain the 
equations 



dU P Til 1 



(2) 



2 



and 

- -Q^-tp v px = (3) 

where t p = (e p — p)/T. For an isotropic system, the 
collision integral I has been evaluated analytically by 
an elaborate approximation 14] . As we cannot assume 
the predominance of normal scattering processes in gen- 
eral cases of anisotropic Fermi liquids of our concern, 
we have to make approximations in another way. To 
simplify the multiple momentum sum in J, we replace 
(l-n(e p > +k ))S(£p+e p >-£ p -k-£p>+k) with (1— n{e p +e p >- 
£p-fe))^(£p'+fe — m)> f° r significant contributions to the col- 
lision term should come from a thermal neighborhood of 
the FS in any case. Moreover, as in the isotropic case, to 



describe the momentum dependence of the solutions 
we decouple the crystal momentum variable p into the 
radial (energy) direction t p = (e p — p)/T and_the per- 
pendicular component J7 p , and set li tP = Mi(Q p )Ni(t p ) 
(i = 1,2). Then the momentum sum is written as 
= T J pdt J dCl p , where p represents the density of 
quasiparticle states (DOS) under the normalization con- 
dition J dftp = 1 . By integrating over the energy vari- 
ables, we obtain equations for Mj(fi p ). 

First we note that we can reproduce the same conduc- 
tivity formula as derived and discussed previously [l3l [l5| . 
By setting Ni(t p ) as a constant, the energy integrals 
which appear in the four terms in I[li] of Eq. @ can 
be calculated analytically. As a result, we obtain 



p 2 T 2 / dn p , / dn p »W^ p „(Mi,p + Mi >p / -M hp „ -Mi, p+I ,<_p<0<$(ep+p'_p» -ft) 
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W p ^ p ,,{Mx^ + Miy - Miy - Mi^ p+ p^pn)p p , p p „ p p+p , _ p „ 



(4) 
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In Eq. ([J]), we used p p — S(e p — p). Assuming the an- 
gular dependence Mi jP oc v px , we obtain the resistivity 
coefficient 

_ Ep.p'.p" W^ p „ PpPp, p p n Pp +p , -p»V p!c {Vp X + V p , x - 



dt t*+ir* J U P X 

which we obtain the thermal resistivity Kq 1 
is remarked that the difference in the integrand kernels 
of Eqs. ([!]) and (O stems from the fact that N 2 (t) for kq 
is an odd function, while N±(t) for a is even. To evalu- 
ate the Lorenz number L = k/uT concretely, assuming 



2,p 



above, 



we obtain 



U p " x 



Vp+p' —p" x) 



4e : 



(Ep Pp V px 



for the electrical resistivity a 1 = AT 2 . 

On the other hand, we have to make a further approx- 
imation for kq. By taking the i-derivative at t p = of 



(8) 



Eq. ©, while adopting N 2 (t) oc t/(t 2 + n 2 )[ 
obtain similar equations, namely, 



we finally 



v px = p 2 T 2 J dn p , J dCl pll WPP p „(M 2 .p-c(M 2 ^ 
+M 2 , P » + M %p+p ,_ p ,,))8{e p+p'—p" 

-A*), (7) 

where 

WtoothM 2 ^r 2) ^0.162, 
2 / sinhi 



where the Boltzmann constant k is written explicitly [27|. 
The constant prefactor in Eq. (8J is separated as it comes 
from the energy integral, and corresponds to the result 
first derived by Herring; 8], though our result differs by 
7r 2 /12. In addition, the dimcnsionless factor C in Eq. ([5]) 
originates from the directional dependence of scattering, 
and is given explicitly by 
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(9) 



(10) 



Note that one would find a trivial result £ = 1, a con- 
stant Lorenz ratio, if one neglects both of the three vertex 
correction terms following v pix in the parentheses of the 
denominator as well as the numerator of Eq. To have 
a nontrivial effect £ ^ 1, it is essential not to disregard 
detailed momentum structure of the relevant quasiparti- 
cle scatterings on the FS. In fact, £ quantifies the effect 
of FS geometry on the availability of phase space for the 
quasiparticle scatterings to relax the transport currents. 
Essentially, the above expression represents the fact that 
thermal and electrical resistivities are determined mainly 
by different types of scattering processes, that is, normal 
processes are important for the thermal resistivity B, 
while the electrical A is caused by Umklapp processes. 
In fact, for the numerator, or the resistivity coefficient 
A, to take a finite value, there must exist at least a set of 
four momenta Pi (i — 1, 2, 3, 4) on the FS satisfying the 
Umklapp condition pi +p 2 = P3 + P4 + G, where G(^ 0) 
is a reciprocal lattice vector. This is not met if the FS is 
too small. Indeed, there is a limit Fermi radius estimated 
by the relation |pi| = |G|/4 which holds at the threshold 
where all the vectors are parallel or antiparallel. In the 
vicinity of the threshold, one should expect a strong FS 
dependence of £. This is a matter of our concern in the 
following. 



III. RESULTS 
A. Simple Model 

The factor £ may be sensitive to the momen- 
tum dependence of the scattering probability W£ p = 
f (lA^J 2 + \\Afl k f) , where < fc and A^ k arc the 

scattering amplitudes for quasiparticles with parallel and 
antiparallel spins. In particular, the momentum depen- 
dence can give rise to a conspicuous effect in the vicinity 
of quantum critical points where there are quantum fluc- 
tuations localized in k space. To take this into account, 
we may assume 
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l + £ 2 7(k-Q) 



(11) 



FIG. 1: For a parabolic band on a square lattice, the dimen- 
sionless Lorenz factor C in Eq. is shown as a function of 
the filling n per orbital. For £ = (solid line), for f = 1 and 
Q = (tt, 7r) (dashed line, AFM), and for £ = 1 and Q = (0, 0) 
(dotted line, FM). In all cases, £ must vanish at the threshold 
n c 1 ~ 0.2. 



in order to represent the effect of fluctuations peaked 
around k = Q. In addition to the correlation length £, 
one may adopt 7(k) oc J^d (l — e lk d ), where the sum is 
taken over the nearest neighboring lattice vectors d, as it 
gives a simple lattice periodic function to give 7(k) ~ k 2 
in the long wavelength limit k — > 0. 

To illustrate how £ varies, let us first investigate a 
simple two-dimensional (2D) model. We should make a 
special remark, however. In a strictly 2D system, one will 
find that normal forward scatterings make the denomi- 
nator B of Eq. ([9]) logarithmically divergent, as so for 
the inverse lifetime of quasiparticle which corresponds to 
the denominator of Eq. © with c = 0[l6j. On the other 
side, the numerator A still remains finite as noted by Fuji- 
moto et aZ.|17|. In effect, the divergence is suppressed by 
a small decay rate T assigned to the quasiparticle states, 
p p = S(e p — fi). In a real system, such a cutoff must 
be provided by an inevitable effect of three dimension- 
ality of the system or by a finite density of impurities. 
As the dependence on T is logarithmic and weak numer- 
ically, here we present a typical behaviour assuming F as 
a given constant for simplicity. 

First we discuss the simplest case of a parabolic band 
on a square lattice. It is easy to show that the re- 
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FIG. 2: £(= A/B) of a parabolic band on a sc lattice and a fee 
lattice. In the inset, the resistivity coefficients A cx <r -1 /T 2 
and B for the sc lattice are shown as a function of n. 



suit is independent of the quasiparticle mass, so that 
we may simply use e p — p 2 . We obtain Fig. [T] for 
r = 0.02, which representatively shows £ as a function 
of the electron number n per orbital for three types of 
W pp according to Eq. (fTTj) . As expected, characteristic 
anomalies are clearly observed. The solid line for £ = 
(W pp =const.) typically shows the first onset of the 
Umklapp processes involving the smallest reciprocal lat- 
tice vector G = (2tt, 0) at n Cj i = 7r/16 ~ 0.20, as well 
as the second one at n Cj 2 — 0.39 for G = (2ir, 2ir). It 
is noted that the threshold fillings can be easily evalu- 
ated as they are geometrically determined by the given 
FS. The dashed line indicates that the latter structure for 
n > n c 2 is particularly emphasized by the commensurate 
antifcrromagnetic (AFM) fluctuations with Q = (tt,it). 
On the other hand, the ferromagnetic (FM) fluctuations 
with Q = (0, 0) generally suppress £, as they strengthen 
the relative weight of the normal processes contributing 
to the thermal resistivity B. 

Similarly, one may obtain results for three dimensional 
systems, in which no cutoff is required. In Fig. [2l we 
observe the lattice-structure dependence of C for the 
parabolic band in a simple cubic (sc) lattice and a face 

centered cubic (fee) lattice with W pp =const. For the 
sc lattice the first threshold lies at n Cj i ~ 0.065 due 
to G = (2tt,0,0), while it is at n cA ' ~ 0.085 with 
G = (2-7T, 2-7T, — 2tt) for the fee lattice. For the latter, the 
secondary kinks expected at n Cj2 ~ 0.13 and n Ci3 ~ 0.37 
corresponding to G = (47r,47r,0) and (47r,0, 0), respec- 
tively, are not so conspicuous as that found clearly at 
n Cj 2 — 0.19 due to G = (2ir, 2ir, 0) for the sc lattice. In 
the inset of Fig. [2j the resistivity coefficients A and B 
are shown respectively for the sc lattice. We remark that 
the presence of the threshold n C: i may be more easily 
anticipated from a relatively gradual n-dependence of C 
than from the electrical resistivity coefficient A, which 
drops abruptly at n c ,i- The results exemplify that the 
ideal Lorenz ratio is not a constant number but shows 
the lattice structure dependence interestingly. 



TABLE I: Calculated values for the ideal Lorenz ratio are 
compared with experiments. 

C L (1Q- 8 V 2 /K 2 ) L cxp (1Q- 8 V 2 /K 2 ) 
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FIG. 3: For a single band tight-binding model of Na I Co02, 
£ as a function of x as well as the Fermi surfaces at x = 0.6 
are shown to indicate the strong effect caused by a slight 
deformation due to the hopping integrals t2 and £3. 



B. Transition metals 



In principle, we can evaluate C for real materials 
though it would generally require a hard task numeri- 
cally. Let us evaluate them for transition metals in the 
same manner as described elsewhere (l8j. i.e., from the 
result of a first principle band calculation, we pick up 
a main band with the largest DOS, for which we apply 
Eq. (fTU|) . Table Q] shows calculated values along with ex- 
perimental results 0i S, IE IE 0- For these to be evaluated, 
we regarded W£ p as a constant, except Pd for which pre- 
sented also is the result with a paramagnon ferromagnetic 
correlation effect with the correlation length £ = 5A[lj|, 
which is taken into account as above in Eq. (fTTj) . Con- 
sidering the approximations made to derive Eq. (|10[) . for 
the scattering amplitudes and so on, we conclude that 
we could explain a small L cxp of tungsten, among others. 
This must be primarily due to a peculiar FS[20]. For 
platinum, however, the observed value[5j is inexplicable 
by our single band result. 
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C. Na^CoOz 

Lastly, we discuss a tight-binding model of Na x Co02. 
Though this material has attracted much interest by its 
large thcrmopower, it is of particular interest for us be- 
cause of its simple cylindrical FS. As a nearly 2D sys- 
tem on a triangular lattice, it has an almost circular hole 
surface centered at the T point with a slight hexagonal 
distortion bul ging into the sides of the hexagonal Bril- 
louin zone[TTI. l2ll. |2^| . Therefore, it provides us with an 
unique case to realize the nontrivial effect as discussed 
above. Indeed the FS suggests that the Umklapp thresh- 
old must exist around x c ~ 0.6, and one should expect 
that even a slight change in the shape of the FS around 
the threshold could have a striking effect on the doping 
dependence of the genuine quasiparticle transport prop- 
erties and the ideal Lorenz factor. In fact, if we assume a 
undistorted parabolic band, it is straightforward to show 
that we should obtain x c — 2n C) x — la 0.55. Never- 
theless, experimentally, the T 2 dependence of the electri- 
cal resistivity has been observed up to x ~ 0.7[H, |24j . 
around which, therefore, it must be necessary to take into 
account a three dimensional lattice distortion [2l|,[22j ■ For 
definiteness and simplicity, we restrict ourselves to the 
doping regime where the three dimensional effect is ir- 
relevant. We do not consider incipient ferromagnetic in- 
plain correlations [1^, as they would not modify the result 
qualitatively. Below, we pay special attention to an ex- 
pected nontrivial behavior that the factor C should tend 
to vanish as we approach a threshold x < x c . We present 
results for a 2D tight-binding model on a triangular lat- 
tice, in which up to the third neighbor hopping integrals 
t{ (i = 1, 2, 3) are considered as in Ref. [26j, where Kuroki 
and Arita have discussed the thermopower and the elec- 
trical conductivity by a relaxation time approximation 
with a single time scale r. Here the approximation must 
be abandoned from the outset, since otherwise one would 
only obtain a trivial result of a constant value regardless 
of the FS, as mentioned above. 

To see how C tends to vanish around x < x c quali- 
tatively, we show Fig. [3] in which the x-dependence of C 
around x < x c is shown for t\ > and T/ti = 0.02, along 
with the Fermi surface for x = 0.6. In comparison with 
the dashed line for — t$ = 0, the solid line shows the 
strong effect around x < x c due to a slight modification 
of the FS caused by t-i and £3. In effect, as the threshold 
value x c ~ 0.6 itself depends sensitively on portions of 
the FS closest to the Brillouin zone boundary, the result 
cannot be regarded as a quantitative prediction. Fur- 
thermore, unfortunately, it is not easy to compare this 
nontrivial prediction directly with experimental results, 
because the intrinsic electron term of the thermal resis- 
tivity (oc BT) for the perovskite oxides has been com- 
pletely outweighed by contributions due to phonons and 
impurities [23L 1241] . Those extrinsic terms have to be sep- 
arated out properly to verify the nontrivial filling depen- 
dence. 

Lastly, as the vanishment of the factor C essentially 




FIG. 4: As in Fig. |3l shown are the filling x dependences of 
the factor F, proportional to the KW ratio A/7 2 . 



reflects that of the electrical resistivity, it would be in- 
teresting also to calculate the x dependence of the resis- 
tivity coefficient A, or the Kadowaki Woods (KW) ratio 
A/7 2 . To evaluate the resistivity coefficients themselves, 
however, we need the absolute value of the scattering am- 
plitude. Assuming the strong coupling p 2 W p v _ v „ a 7r[l8l|. 
e.g., then we obtain A/7 2 = 9F/167re 2 with 

PpiPP2PP3PP3-Pl-P2 V Plx( V Plx + V Plm ~ V P3IC ~ V Pl+P2~P3,x) 



1,2,3 




F = 



As in Fig. [3] for C, we show F as a function of x in 
Fig. 21 The ratio A/7 2 also vanishes as x approaches 
expected. However, numerically, we find that it 
should remain of the order of a common value, A/7 2 ~ 
1 x 10~ 5 ,uf2 cm(mol K/mJ) 2 . Li et al. [23| have observed 
a strongly enhanced deviation from this standard value 
for x = 0.7, which should be well beyond the threshold 
x c . It is difficult to regard the numerical deviation simply 
as a FS effect. Apart from taking into account a relevant 
three dimensional effect, one would have to assume an 
enhanced scattering p 2 W^ v _ pll 3> 1, which may be due to 
strong scatterings caused by a proximity to some sort of 
instability [H. 



IV. SUMMARY 

In summary, on the basis of anisotropic Fermi liquid 
theory, we investigated the ideal Lorenz ratio for cor- 
related metals by taking due care of the momentum de- 
pendence of transport relaxation processes due to mutual 
elastic scatterings between quasiparticles. It was shown 
explicitly that the ideal Lorenz ratio of a correlated elec- 
tron system is not a constant, but may vary drastically 
in circumstances and even be made vanishingly small, 
not only by quantum fluctuations but by a filling control 
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around thresholds of Umklapp scattering channels. Al- 
though it might not be easy to extract the ideal Lorenz 
value in practice, theoretically we pointed out and dis- 
cussed that such a nontrivial effect should be expected 



in such a simple single-band system as Na^CoC^- 

The numerical calculations were carried out on Al- 
tix4700 at Shizuoka University Information Processing 
Center. 
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